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Abstract. The spectrum of thermal gravitational waves is obtained by including the 
high frequency thermal gravitons created from extra-dimensional effect and is a new 
feature of the spectrum. The amplitude and spectral energy density of gravitational 
waves in thermal vacuum state are found enhanced. The amplitude of the waves get 
O . modified in the frequency range (10~ 16 -10 8 Hz) but the corresponding spectral energy 

density is less than the upper bound of various estimated results. With the addition 
£h ■ of higher frequency thermal waves, the obtained spectral energy density of the wave in 

thermal vacuum state does not exceed the upper bound put by nucleosynthesis rate. 
The existence of cosmologically originated thermal gravitational waves due to extra 
dimension is not ruled out. 



PACS numbers: 98.70.Vc,98.80.cq,04.30.-w 
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1. Introduction 

One of the predictions of general theory of relativity is the existence of gravitational 
waves. The sources of generation of these waves vary from the dynamics of early universe 
to massive astrophysical objects such as neutron star binaries and black hole mergers 
etc. Thus the waves have a wide spectrum of frequencies vary from very low to high 
(9(10~ 19 Hz -10 10 Hz). It is possible to discriminates the relic waves from other sources 
on the observational point of view also. The relic gravitational waves are paramount 
importance in cosmology because it provides valuable information on the conditions of 
very early universe. It is believed that the relic gravitational waves are mainly generated 
during inflationary epoch. And the waves that amplified during the inflation are low 
frequency only. Since the higher frequency waves are outside the "barrier" (horizon) [TJ 
the corresponding amplitudes decreased during the evolution of the universe. The 
features of relic gravitational waves of very high frequency range is interesting though 
the energy scale of conventional inflationary models are not favoring for it. However 
if extra dimensions exist (for a review, see [2] and motivation for extra dimensions 
[3]) the graviton background can have a thermal spectrum [4]. According to the 
extra dimensional models the thermal gravitons with very high frequency range also 
be observed with a specific peak temperature today jl] and therefore the detection 
of very high frequency thermal gravitational waves is an interesting test to see the 
possibility of existence of extra dimensions as well. These thermal gravitational waves 
can contribute to the higher frequency range of the spectrum. The existence of thermal 
graviton background with the black body type spectrum is also discussed in [6]. If 
the inflation was preceded by a radiation era, then there would be thermal gravitational 
waves at the time of inflation [6] . The generation of tensor perturbations during inflation 
by the stimulated emission process leads to the existence of thermal gravitational waves 
[7]. Direct detection of the thermal gravitons is challenging but may be possible in the 
near future with the 21-cm emission line of atomic hydrogen. 

The inflationary scenario [8] predicts a stochastic cosmic background of 
gravitational waves (CGWB) |5J. The spectrum of these relic gravitational waves 
depends not only on the details of expansion during the inflationary era but also the 
subsequent stages, including the current epoch of the universe. Computation of the 
spectrum of the waves for matter dominated universe is usually done in decelerated 
expanding model [9]-[Tl]. The resulting spectrum is used for putting constrain on the 
detection of gravitational waves originated from sources other than early universe epoch. 
The result of astronomical observations on SN la [15]- [16] shows that the universe is 
currently under going accelerated expansion indicating a non-zero cosmological constant. 
According to the ACDM concordance model, the observed acceleration of the present 
universe is supposed to be driven by the dark energy. Effect of the current acceleration 
on the nature of the spectrum and spectral energy density of the relic gravitational waves 
is studied (I0],[l7]. And shown that the current acceleration phase of the universe does 

X the terminology "barrier" is adopted from pQ. 
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change the shape, amplitude and spectrum of the waves [T7] . 

In the present work, we consider contribution of very higher frequency relic thermal 
gravitational waves to its spectrum and spectral energy density for the decelerated as 
well as accelerated universe. The focus of the present work is on the spectrum of the 
higher frequency range of the waves due to extra dimensional effects. The normalization 
of the spectrum is being done with the measured CMB anisotropy spectrum of the 
WMAP. The inclusion of the higher frequency relic thermal gravitational waves leads to 
enhancement of the spectrum. This enhancement leads to modification of the amplitude 
of the spectrum in the frequency range (10~ 16 -10 8 Hz) as an additional feature and is 
possible to compare these with the sensitivity of Advanced. LIGO (Adv.LIGO), Einstein 
Telescope (ET) and LISA missions. The corresponding spectral energy density can be 
compared with estimated upper bound of various studies. Also can check whether the 
inclusion of higher frequency thermal gravitational waves in the total spectral energy 
density exceed the upper bound of primordial nucleosynthesis rate or not. In the present 
work, we use the unit c = h — ks — 1- 



2. Gravitational waves spectrum in expanding universe 

The perturbed metric for a homogeneous isotropic flat Friedmann-Robertson- Walker 
(FRW) universe can be written as 

ds 2 = S 2 {r]){d V 2 - (Sij + hijjd^dx 3 ), (1) 

where S(rj) is the cosmological scale factor, 77 is the conformal time and Sij is the 
Kronecker delta symbol. The are metric perturbations field contain only the pure 
gravitational waves and is transverse-traceless i.e; Vj/i*- 7 = 0, S^hij = 0. 

The present study mainly deals with amplitude and spectral energy density of the 
relic gravitational waves generated by the expanding spacetime background. Thus the 
perturbed matter source is therefore not taken into account. The gravitational waves 
are described with the linearized field equation given by 

VcjV^V^x,!,))^. (2) 

The tensor perturbations have two independent physical degrees of freedom and are 
denotes as h + and h x , called polarization modes. To compute the spectrum of 
gravitational waves h(x, rf) in the thermal states, we express h + and h x in terms of 
the creation (a') and annihilation (a) operators, 

' 'alhl( V )^ + at P h^( V )e-^], (3) 



2k 

where k is the comoving wave number, k — |k|, l p i — yG is the Planck's length and 
p = +, x are polarization modes. The polarization tensor e?(k) is symmetric and 
transverse-traceless fc J e?-(k) = 0,5 y e?-(k) = and satisfy the conditions e yp (k)e?- (k) = 
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2<5 pp / and e?-(— k) = e?-(k), the creation and annihilation operators satisfy [flk' a k' p/ ] = 
^pp'^ 3 (k — k'), the initial vacuum state is defined as 

o£|0> = 0, (4) 

for each k and p. The energy density of the gravitational waves in vacuum state is 

pi 

For a fixed wave number k and a fixed polarization state p the linearized wave 
equation (j2j) gives 

K + 2jh' k + k 2 h k = 0, (5) 

where prime means derivative with respect to the conformal time. Since the polarization 
states are same, we here onwards denote hk(rj) without the polarization index. 

Next, we rescale the filed hk(rj) by taking h k {rf) = fk(v)/^( r ))i where the mode 
functions fk{v) obey the minimally coupled Klein-Gordon equation 

f k '+(k 2 -^)f k = V. (6) 

The general solution of the above equation is a linear combination of the Hankel function 
with a generic power-law for the scale factor S = r] q given by 

f k (v) = A^H^kv) + Bk^hqH^ikri). (7) 

For a given model of the expansion of universe, consisting of a sequence of successive 
scale factor with different q, we can obtain an exact solution fk(v) by matching its value 
and derivative at the joining points. 

The approximate computation of the spectrum of gravitational waves is usually 
performed in two limiting cases depending up on the waves that are within or outside 
of the barrier. For the gravitational waves outside barrier (k 2 ^> S"/S, short wave 
approximation) the corresponding amplitude decrease as hk oc 1/S(rj) while for the 
waves inside the barrier (k 2 <C S"/S, long wave approximation), hk = Ck simply a 
constant. Thus these results can be used to estimate the spectrum for the present 
epoch of universe. 

The history of overall expansion of the universe can be modeled as following 
sequence of successive epochs of power-law expansion. 
The initial stage (inflationary) 

S( V ) = l \v\ 1+ ", -oo<V<Vi, (8) 

where l + /3<0,^<0 and l is a constant. 
The z- stage 

S(n) = S z {r]- Vp ) 1+ ^, r]i<V< Vs, (9) 

where /3 S + 1 > 0. Towards the end of inflation, during the reheating, the equation of 
state of energy in the universe can be quite complicated and is rather model-dependent 
[T8] . Hence this z-stage is introduced to allow a general reheating epoch, see for details 

□n. 
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The radiation-dominated stage 

s (v) = S e (v-Ve), Vs<V<V2, (10) 

The matter-dominated stage 

S(r)) = S m (r) - r/ m ) 2 , r] 2 <r]<r] E , (11) 

where T]e is the time when the dark energy density pa is equal to the matter energy 
density p m . Before the discovery of accelerating expansion of the universe, the 
current expansion is used to take as decelerating one because of the matter-dominated 
stage. Thus, following the matter-dominated stage, it reasonable to add an epoch of 
accelerating stage, which is probably driven by either the cosmological constant, or the 
quintessence, or some other kind of condensate [!§] . The value of redshift ze at t]e 
is (1 + ze) = S(t] )/ S(r] E ), where r] is the present time. Since pa is constant and 
p m {rj) oc S~ 3 (i]), we get 

pA = 2± = i. (12) 

Pmij]E) Pm(^o)(l + ^i?) 3 

If the current value of £7 a ~ 0.7 and fl m ~ 0.3, then it follows that 

i + z E =(- F r) ~i- 33 - (1 3 ) 



The accelerating stage (up to the present) 

s{n) = to\ , n-Va\~ 1 , ve^v^vo- (14) 

This stage describes the accelerating expansion of the universe. And is a new feature 
and hence its influence on the spectrum of relic gravitational waves is of interesting to 
study. It is be noted that the actual scale factor function S(r]) differs from equation (fT4"|) . 
since the matter component exists in the current universe. However, the dark energy is 
dominant, therefore (1141) is an approximation to the current expansion behaviour. 

Given S(rj) for the various epochs, the derivative 5" = dS/drj and ratio S'/S follow 
immediately. Except for (3 S which is imposed upon as the model parameter, there are 
ten constants in the expressions of S(rj). By the continuity conditions of S(r]) and 
S'(r]) at four given joining points r)i,r) s ,r)2, and t/e, one can fix only eight constants. 
The remaining two constants can be fixed by the overall normalization of S and the 
observed Hubble constant as the expansion rate. Specifically, we put \r] — r] a \ = 1 for 
the normalization of S, which fixes the T] a , and the constant £q is fixed by the following 
calculation, 

H~\3' 



(15) 

vo 

where £o is the Hubble radius at present. 

In the expanding Friedmann- Robertson- Walker spacetime the physical wavelength 

Co 



is related to the comoving wave number as A = 2nS ^ ; and the wave number fco 



corresponding to the present Hubble radius is ko = 2itS ^ = 2n. And there is another 
wave number fog = 2?r ^ E ^ = , whose corresponding wavelength at the time t\e is 
the Hubble radius 1/H. 
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By matching S and S' / S at the joint points, one gets 

/ = « s (2+/3) c 2 v cfcr l? , (i6) 

where b = |1 + f3\~^\ which is defined differently from [20], Ce = f§§, C2 = 
C s = §7^4, and Ci = ffej- With these specifications, the functions S^) and S'(r))/S(rj) 
are fully determined. In particular, S'(t))/S(t)) rises up during the accelerating stage, 
instead of decreasing as in the matter-dominated stage. This causes the modifications 
to the spectrum of relic gravitational waves. 

3. Gravitational waves spectrum in thermal vacuum state 

The power spectrum of gravitational waves is defined as 

rlk 

/ i 2 (A;,r / )^ = (0|^(x,r ? )^(x,r ? )|0), (17) 







Substituting equation (J3]) in ffTTj) and taking the contribution from each polarization is 

same, we get 

41 , 

h(k,r)) = ^k\h(r])\ . (18) 



Thus once the mode function h{rj) is known, the spectrum h(k, if) follows. 

The spectrum at the present time h(k, rj ) can be obtained, provided the initial 
spectrum is specified. The initial condition is taken to be the during the inflationary 
stage. Thus the initial amplitude of the spectrum is given by 



k 



2+j3 



h(k,T H )=A[- j -\ , (19) 

where A = 8y/7r 1 -^- is a constant. The power spectrum for the primordial perturbation 
of energy density is P(k) oc | h(k,f]o) | 2 and in terms of initial spectral index n, it is 
defined as P(k) oc k n ~ l . Thus the scale invariant spectral index n — 1 for the pure de 
Sitter expansion can be obtained with the relation n = 2/3 + 5 for 0— - 2. 

An effective approach to deals with the thermal vacuum state is the thermo-field 
dynamics (TFD) [21] . In this approach a tilde space is needed besides the usual Hilbert 
space, and the direct product space is made up of the these two spaces. Every operator 
and state in the Hilbert space has the corresponding counter part in the tilde space [2Tj . 
Therefore a thermal vacuum state (Tv) can be defined as 

\Tv) = T(e k )\0 0), (20) 

where 

T{9 k ) = exp[-# fc (a k a k - a^a^)], (21) 

is the thermal operator and |0 0) is the two mode vacuum state at zero temperature. 
The quantity 6 k is related to the average number of the thermal particle, n k = sinh 2 ^^. 
The fik for given temperature T is provided by the Bose-Einstein distribution = 
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[exp(k/T) — l] -1 , where Uk is the resonance frequency of the field. The a^, a k and Ok, 
a/ k , are respectively the annihilation and creation operators in Hilbert and tilde space, 
satisfy the usual commutation relations, [a^, a h ,] = [5k, a^,] = 5 3 (k — k') . And all other 
commutation relations of these operators are zero. By the appropriate action of the 
operator (1211) on and a k , we get [22] 



(22) 



7~WT = cik cosh 9k + a k sinh 9k, 



-Pair 



it cosh 9k + ak sinh 9 



Hence the occupation number in thermal vacuum state can be written as 



(4 a k') 



i 



,fc/T _ 1 



5 3 (k-k') 



(23) 



Thus, using Eq.flBj) and Eqs. (12011231) in Eq.ffTTl) the power spectrum in thermal vacuum 
state is obtained as 



hi 



(k,T)) 



16L 2 



7T 



-k 2 \ h(r)) | 2 coth 



A 

2T 



Thus in comparison with Eq. ffl9l) . the spectrum is expressed as 



h(k,r}i 



A 



k 



2+/3 



f ] coth 1 / 2 
k 



2T 



(24) 



(25) 



The last term becomes significant when the ratio k/(2T) is less than unity. The 
wave number k and temperature T are comoving quantity which are related to the 
physical parameters at the time of inflation, see for details [6]. Thus it is expected an 
enhancement of the spectrum by a factor coth 1,/2 [/c/2T] = coth 1//2 [H Si/2Ti\ . 

It is convenient to consider the amplitude of waves in different range of wave 
numbers [T7]. Thus the amplitude of the spectrum in thermal vacuum state for different 
ranges are given by 

(i) when k < k E , the corresponding wavelength is greater the present Hubble radius. 
Thus the amplitude remain as the initial one and can be written as 



h T {k,r} ) = A 



k\i+P 



coth 1/2 



L2T 



(26) 



(ii) the amplitude remains approximately same as long as the wave inside the barrier 
but begin to decrease when it leaves the barrier by a factor 1/S(rj), depending the value 
of scale factor at that time. This process continue until the barrier becomes higher than 
A; at a time rj earlier than 770, so the amplitude has decreased by the ratio of the scale 
factor at the time of leaving the barrier 5& to its value at 77, S(j]). This is in the range 
k E <k < ko- 



h T {k,T)o) = A 



k 



coth 1/2 



L2T 



ze) 



(27) 



Note that this range is a new feature on account of the current acceleration of the 
universe which is absent in the decelerating model as pointed out in [17] . The amplitude 
of the waves that left the barrier at S& with waves numbers k > k has been decreased 
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up to the present time by a factor Sb/S(rj ). This affect the amplitude of the present 
spectrum and is obtained as 



k \ 2 +P 



k i S, 



2T 



h T (k, Vo ) = A(-) coth 1 / 2 — — \. (28) 



S(vo) ' 



This result can be used to obtain the spectrum of the waves in the remaining range of 
wave numbers. 

(iii) the wave number that does not hit the barrier in the range k < k < k 2 gives 
the amplitude as follows 

MM.H^)W/>[A]_J_, (29) 

the spectrum in this interval is differ from that of the matter dominated case by a the 
factor r 1+ZE \3 ■ The wave lengths of the spectrum in the range are long but smaller than 
the present Hubble radius. 

(iv) in the range of wave number k 2 < k < k s , the amplitude is 

This is the interesting range on the observational point of view of Adv.LIGO, ET and 
LISA. Note that the temperature dependent factor in this range is negligible hence the 
term is dropped out because of the low temperature nature of the relic waves. 

(v) for the wave number range k s < k < k\ which is in the high frequency case and 
gives the corresponding amplitude as 

h T (k, m )=A{-) (_) (^)cothV L2T 



12 



:i+^) 3 ' 

(31) 



In the usual case the temperature dependent term can also be neglected however 
the extra dimensional scenario predicts higher temperature for the thermal gravitational 
waves, hence the term again becomes significant. Therefore the contribution from 
the thermal relic gravitational waves is expected increase the amplitude of spectrum 
particularly in the higher frequency range also. 

It is to be noted that in (iv) the thermal contribution in k 2 < k < k s range 
is negligible due to the temperature dependent term. Similarly the thermal effect is 
insignificant in the range k s < k < k\ also. However by taking into account the 
extra dimensional effect, the spectrum of relic waves is peaked with a temperature 
T*=1.19 x 10 25 Mpc -1 [3] (See, appendix A for a brief discussion on T* from extra 
dimensional scenario.). Therefore it is expected enhancement for the amplitude of 
spectrum (orange lines, Figs.jl] and j4]) in the range k s < k < k\ compared to T = 
case for the accelerated as well as decelerated universe. But at the same time, ignoring 
the thermal contribution to the amplitude of spectrum in the range k 2 < k < k s leads 
to a discontinuity at k s , see Fig. [I]. This is evaded by fitting a new line in the range 
k 2 < k < k s for the amplitude h of Eq.( |30l) as follows. 
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Let the amplitude of the wave in the range k < k < k 2 is given by 
be rewritten as 



h 1T (k,T]o) = A 



k\P 



coth 1 / 2 



L2TJ (l + z E ) 



and can 



(32) 



and the amplitude in the k s < k < k\ is given by f l3Tj) also rewritten as 



h 2T {k,r] ) = A 



ko 
k 2 ) 



Wh 1/2 



2Tj (1 + z e ) 



(33) 



k ) \k , 

Thus the new slope for Eq.f l30|) . in the range k 2 < k < k s , can be obtained by taking 
y = log 10 (/i) and x = log 10 (fc), then 

log M W - logioWi = ^^T^m ^sM - logio(^)), (34) 

where the subscribes z and / are respestively indicating the first and last points of the 
straight line. By putting ki = ^2 from Eq. (l32|) and kf = k s from Eq.( )33|) in Eq. (l34p . we 

get y 

h = (h 1T ) k2 g(k), 



where 



and 



g(k) 



k \7 
k 2 , 



7 



logi (^2T)fc s - logi (/lll\ 



fe 2 



lo gio(^) -log 10 (fc 2 



lo Sio((|) 


1 /9 

coth 1/2 


' k s " 
2T* 


) 


lo 


Sio \T 2i 





(35) 
(36) 

(37) 



is the slope of the line and thus we find the amplitude, for convenience we call it as 
'modified amplitude', given by 



h(k,rj ) = A 



ko\P 



1 / k \7 

'i+z E )A~k 2 , 



(3f 



When becomes zero Eq. fl37|) leads to 7 = 1 + (3, and hence ( 1301 is recovered from 
f )38|) in the range k 2 < k < k s . 

The overall multiplication factor A in all the spectra is determined in absence of the 
temperature dependent term with the CMB data of WMAP [17] . This is based on the 
assumption that the contribution from gravitational waves and the density perturbations 
are the same order of magnitude or if the CMB anisotropies at low multipole are induced 
by the gravitational waves, therefore it is possible to write AT/T ~ h(k,T]o)- The 
observed CMB anisotropies [23] at lower multipoles is AT/T ~ 0.37 x at I ~ 2 
which corresponds to the largest scale anisotropies that have observed so far. Thus 
taking this to be the perturbations at the Hubble radius gives 

h(k , no) = A 1 = 0.37 x KT 5 . (39) 

(l + z E y 



1 /9 

here, coth ' 



2T 



= 1. 
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Figure 1. The amplitude of the gravitational waves for the accelerated (solid lines) 
and decelerated (dashed lines) universe. 

However, there is a subtlety in the interpretation of AT/T at low multipoles, whose 
corresponding scale is very large ~ t§. At present the Hubble radius is £q, and the Hubble 
diameter is 2£ . On the other hand, the smallest characteristic wave number is &e, whose 
corresponding physical wave length at present is 2irS(rjo)/kE = ^o(l + Ze) — l-32£ , 
which is within the Hubble diameter 2£ , and is theoretically observable. So, instead of 
Eq. (l39~l) . if AT/T ~ 0.37 x 10~ 5 at I ~ 2 were taken as the amplitude of the spectrum 
at ve, one would have /it(&e,?7o) = + ze) 2+i3 = 0.37 x 10~ 5 , yielding a smaller A 
than that in Eq. fl39|) by a factor (1 + zE) 1 ^ 13 ~ 2.3 [17J. The allowed values of (3 and 
(3 S are obtained and are respectively give by (5 = —1.9, and (3 S = —0.552 [T7] . 

Next, we obtain the spectrum in the thermal vacuum state with the following 
parameters. By taking k = 2nv, v E = 1.5 x 10 -18 Hz, u = 2x 10 -18 Hz, u 2 = 117x 10~ 18 
Hz, u s = 10 8 Hz, Ui = 3 x 10 10 Hz ( the value of v\ is taken such a way that spectral 
energy density does not exceed the level of 10~ 6 , as required by the rate of primordial 
nucleosynthesis). The range of frequency is chosen in accordance with generation of 
gravitational waves that vary from early universe to various astrophysical sources. And 
hence the range is matching with the interest of CMB, Adv.LIGO, ET and LISA 
operations for detection of the gravitational waves. The spectrum is computed in the 
thermal vacuum state with the chosen values of the parameters for the accelerated as well 
as decelerated model with T = O.OOlMpc -1 in the low range k < k 2 - This temperature 
is considered in the context of B mode of CMB spectrum in thermal state [6]. And T* 
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log io v(Hz) 

Figure 2. Comparison of the modified amplitude of the spectrum for the accelerated 
(solid black and green lines) and decelerated (dashed black and green lines) universe 
with the sensitivity curves of Adv.LIGO [24] and ET [25] . 

= 1.19x 10 25 Mpc _1 [Jj] for the high range k s ^ k ^ k\ which is from the extra dimensional 
scenario [4]. Since we use the natural unit, the wave number and temperature that 
appear in the temperature dependent term of the spectrum is computed numerically in 
the Mpc" 1 unit. The obtained spectra are normalized of the CMB anisotropy spectrum 
of WMAP data. The amplitude of the spectrum of the thermal gravitational waves 
is enhanced compared to its zero temperature case (vacuum case). It is observed that 
the spectrum for T = O.OOlMpc -1 get maximum enhancement ~ 1.51 times than the 
vacuum case, at k — k^^ and it is ~ 20 times for T i(i =1.19xl0 25 Mpc~ 1 at k = k s . 

The plots for the amplitude of spectrum hT(k,r] Q ) versus the frequency v for 
= —1.9 and (3 S = —0.552 are given in Fig. [I]. The amplitude of the spectrum get 
enhanced in the frequency ranges, 10 -19 Hz< v < 1.49 x 10 _17 Hz, and v s < v < v\ 
( the lower value of this range is selected such way that thermal enhancement of the 
spectral density does not exceed the upper bound of nucleosynthesis rate.) due to the 
thermal effect of gravitational waves but for the range 1.49 x 10 _17 Hz < v < v s there is 
a suppression because of the coth 1 ^ 2 [/c/2T] term. For comparison, the amplitude of the 
spectra are plotted for the decelerated and accelerated universe, see Fig. p] . 

The new enhancement of the gravitational waves spectrum due to the extra 
dimensional effect (the modified amplitude, see Fig. [TJ, light green lines) can be 
compared with the sensitivity of Adv.LIGO, ET and LISA. An analytical expressions 

|| here, T* = 0.905 K= LMxlO^Mpc- 1 . 
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Figure 3. Comparison of the modified amplitude of the spectrum for the accelerated 
(solid black and green lines) and decelerated (dashed black and green lines) universe 
with the LISA sensitivity curve. 



for the Adv.LIGO and ET interferometers are discussed in [26]. For Adv.LIGO and ET 
cases, consider the root mean square amplitude per root Hz which equal to 

M. (40) 
sju 

The comparison of the sensitivity (10 Hz - 10 4 Hz) curve of the ground based 
interferometer Adv.LIGO [24] with the gravitational wave spectra of (5 = —1.9 for the 
accelerated and decelerated universe are given in Fig. [2]. Thus it shows that the Adv. 
LIGO is unlikely to detect the enhancement of the spectrum from the extra dimensional 
effect with its current stands but be possible with the sensitivity of ET. 

Next, we compare the enhancement of the spectrum with the sensitivity (10 -7 
Hz - 10° Hz) of space based detector LISA [27J. It is assumed that LISA has one 
year observation time which corresponds to frequency bin Au = 3 x 10 ~ 8 Hz ( one 
cycle/year) around each frequency. Hence to make a comparison with the sensitivity 
curve, a rescaling of the spectrum h(u) is required in Eq.l fTSl) into the root mean square 
spectrum h{y, Au) in the band Au, given by 



h(u, Au) = h(u)J^-. (41) 

The plots of the LISA sensitivity with the modified amplitude of the spectrum are given 
in Fig. [3] for the accelerated and decelerated universe. This show that the LISA is 
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unlikely to detect the spectrum with the new enhancement feature of the gravitational 
waves. 

The spectral energy density parameter £l g (u) of gravitational waves is defined 
through the relation p g /p c = f Q g (u)^-, where p g is the energy density of the 
gravitational waves and p c is the critical energy density Thus 

n g (v) = ^hUk,Vo)(^) 2 - (42) 

Since the spacetime is assumed to be spatially flat K — with Q — 1, the fraction density 
of relic gravitational waves must be less than unity, p g /p c < 1- After normalization of 
the spectrum by using Eq.(|39l). we integrate J Q, g {y)dvjv from = 10~ 19 Hz up to 
v\ = 3 x 10 10 Hz, with j3 = —1.9 and (3 S = —0.552. The integral is evaluated for the 
thermal case and zero temperature case, the obtained results for the accelerated universe 
are 

(a) is* < v < v E , 

0. 

0.001 Mpc-\ 
0. 

0.001 Mpc~\ 
0. 

0.001 Mpc-\ 
0. 

(e) u s < v < z/i, 

^ = 2.7xl0~ 6 , T = 0. 

Pc 

^ = 6.67 x 10~ 6 , T, = 1.19 x 10 25 Mpc' 1 . 

Pc 

It is to be noted that in (d) the thermal case are not shown because the thermal 
contribution in this frequency range is negligible due to the temperature dependent term. 
However by taking into the extra dimensional effect, the upper limit of temperature of 
the relic waves is to be T*=1.19x 10 25 Mpc _1 . Thus it is expected an enhancement of 
the spectral energy density in range v s < v < V\ compared to T = case for the 
accelerated as well as decelerated universe. But at the same time ignoring the thermal 



^ = 5.8xl0- n , T = 

Pc 

^ = 8.8xl0" n , T = 

Pc 

(b) v E < v < v H , 

^ = 2.3xl0" n , T = 

Pc 

^ = 3.5xl0- n , T = 

Pc 

(c) v H < v < V2, 

^ = 2.4xl0- n , T — 

Pc 

^ = 3.7xl0" n , T = 

Pc 

(d) u 2 < v < u s , 

^ = 8.97 x 10~ 9 , T = 
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Figure 4. The spectral energy density of the gravitational waves for the accelerated 
(solid lines) and decelerated (dashed lines) universe. 



Table 1. Comparison of the estimated upper bound of spectral energy density 
of various studies with the present work. Here rig dec ' and Vlg^ are respectively 
the spectral energy density of the relic gravitational waves in the decelerated and 
accelerated universe of the present study and n g est ^ is the estimated upper bound of 
various studies. 



Frequency(^) Hz 


n ( g dec \u) 


tt g acc \v) 




10" 9 - 1(T 7 


4.98 x 1CT 9 


1.03 x 10" 9 


2 x 10" 8 [25] 


69 - 156 


34.84 x 10" 8 


7.2 x 10" 8 


8.4 x lO" 4 ^ 


41.5- 169.25 


4.93 x 1CT 7 


1.02 x 10~ 7 


6.9 x 10- 6 [H] 



contribution on the spectral density in the range v 2 < v < v s leads to a discontinuity at 
u s , see Fig. j3]. This problem is solved by fitting a new line as discussed in the context of 
estimation of the amplitude of the spectrum and hence recomputed the spectral density 
in the range z/ 2 < v < v s which gives the new value — = 8.21 xl0~ 7 . This changes the 
slope indicating enhancement of the spectral energy density of the gravitational waves 
in the range v<i < v < u s , green lines, Fig. j3]. 

The enhancement spectral energy density £l g (v) in (d) can be compared with 
the estimated upper bound of various studies and are given in Tab.[T]. Thus Q^g ec?l 
and Qj (v) are less than the upper bound of the estimated values of the respective 
frequency range. 



Modified Amplitude of Gravitational Waves Spectrum 



15 



Further see that the contribution to p g j p c from the low frequency range is 
O(l0~ u — 10 -10 ) while from the higher frequency range it is O(10~ 6 ). Since the order 
of contribution to the total p g j p c from the lower frequency side is very small in contrast 
with higher frequency side, we get for the accelerated universe as 

p g /p c ~ 6.67 x 1(T 6 < v < v u (43) 

and is the same order as that of the zero temperature case. However p g /p c of the 
gravitational waves with T / is higher than the zero temperature case at lower 
frequency range < v < v-i. Therefore it is expected an enhancement for the spectral 
energy density in the thermal vacuum state in the frequency range < v < v 2 
and actually it is the range of interest on the observational point of view of the relic 
gravitational waves. The total estimated value of p g / p c by including the thermal relic 
gravitational waves in the very high frequency does not alter the upper bound of the 
nucleosynthesis rate. Thus the relic thermal gravitons with very high frequency range 
are not ruled out and is testable with the upcoming data of various missions for detecting 
gravitational waves. 

4. Discussion and conclusion 

Gravitational waves are one of the classical predictions of Einstein's general theory of 
relativity. The gravitational waves are generated during the very early evolution stages 
of the universe as well as from the various astrophysical objects. Therefore frequency of 
the waves are varying very widely. There are many on going experiments to detect these 
waves and the interested range of frequency is from 10 -19 Hz to 10 10 Hz. The existence 
of gravitational waves with very high frequency range is not favoring by the energy 
scale of the conventional inflationary scenario. However the very high frequency range 
gravitational waves are interesting candidates in the models with extra dimensions. The 
extra dimensional theories predict the existence of thermal gravitons with black body 
type spectrum. These relic thermal gravitational waves can also add to the spectrum of 
the waves thus the corresponding amplitude also get enhanced. The nature of spectrum 
of the waves to be observed today is dependents on the evolution history of the universe. 
Before the result of SN la observations, the current evolution of the universe is used 
to consider as matter dominated with decelerated expansion. But, according to the 
ACDM concordance model the present universe is supposed to be driven by dark energy 
resulting an accelerated phase. If this is the case then the spectral property of the waves 
to be studied by taking into account of the current acceleration of the universe. In the 
present work we mainly considered the very high frequency range (The low frequency 
range thermal gravitational case is carried out by us without including the very high 
frequency thermal waves that comes from extra dimensional scenario and the work 
is under preparation. The enhancement of the lower frequency range is shown with 
red lines, see Figs. fTH4] ) of relic gravitational waves in the thermal vacuum state and 
obtained the spectrum for the accelerated as well as decelerated models. The obtained 
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spectra are normalized with the WMAP data. It is observed that the inclusion of the 
very high relic thermal gravitational waves leads to a discontinuity in the amplitude of 
the spectrum at v s (see Fig.[Tj). This is due to the fact that the temperature dependent 
term is insignificant in the higher frequency side of the range v 2 < v < v s . To evade 
this problem a new equation of line is derived and thus the amplitude get enhanced in 
the range z/ 2 < v < v s . This is the new feature of the spectrum and we designates it 
as the 'modified amplitude' of the spectrum. The modified amplitude of the spectrum 
can be compared with the sensitivity of the Adv.LIGO, ET and LISA missions. The 
comparison of the Adv.LIGO sensitivity shows that the modified amplitude is unlikely to 
detect with its current stands of LISA or the improved sensitivity of Adv.LIGO. Where 
as the proposed sensitivity of the ET is promising to verify the modified amplitude with 
its upcoming mission data. 

The spectral energy density of the gravitational waves is estimated in thermal 
vacuum state for the accelerated and decelerated universe. It is observed that the 
spectral energy density get enhanced in the lower frequency range 0(1O -11 — 10~ 10 ) and 
from the higher frequency range it is (9(10 -6 ). A comparison of the estimated upper 
bound of spectral energy density of various studies with the present work is done. It 
shows that and Vfg CC ^ are less than the estimated upper bound of various studies. 

The total estimated value of p g /p c by including the very high frequency thermal relic 
gravitational waves does not alter the upper bound of the nucleosynthesis rate. Thus the 
relic thermal gravitons with very high frequency range are not ruled out and is testable 
with the upcoming data of various missions for detecting gravitational waves. 



Appendix A 

Extra dimensional Scenario and Thermal Gravitons 



Cosmology with extra dimensions have been motivated since Kaluza and Klein (KK) 
showed that classical electromagnetism and general relativity could be combined in a 
five- dimensional framework. The modern scenarios involving extra dimensions are being 
explored in particle physics, with most models possessing either a large volume or a large 
curvature. Although there exist different models of extra dimensions, there are some 
general features and signals common to all of them. 

In presence of D extra spatial dimensions, the 3+D+l- dimensional action for 
gravity for can be written as 

Rd 



s = J d 4 x J dD y^~9D l ^ r> GD + V-gL r 



where 



G £) — G 



m 



N~ 



pi 



777 



2+D ' 
D 



(44) 



(45) 



and g is the four dimensional metric, G^ is Newton's constant, go,GN and Rd denote 
the higher dimensional counter parts of the metric, Newton's constant, and the Ricci 
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scalar, respectively, nip is the fundamental scale of the extra dimensional theory. 

Since the gravitational interactions are not strong enough to produce a thermal 
gravitons at temperatures below the Planck scale (m p i ~ 1.22 x 10 19 GeV), the 
standard inflationary cosmology predicted the existence the cosmic gravitational waves 
background which are non-thermal in nature. However if the universe contains extra 
dimensions that can generate the thermal gravitational waves and its shape and 
amplitude of the CGWB may change significantly. This can happen when energies in the 
universe are higher than the fundamental scale m^, the gravitational coupling strength 
increases significantly, as the gravitational field spreads out into the full spatial volume. 
Instead of freezing out at ~ 0(m p i), as in 3+1 dimensions, gravitational interactions 
freeze-out at ~ 0(m£>). If the gravitational interactions become strong at an energy 
scale below the reheat temperature (m^ < Trh), gravitons get the opportunity to 
thermalize, creating a thermal CGWB. The qualitative result, the creation of a thermal 
CGWB if m D < T RH , is unchanged by the type of extra dimensions chosen [3]. 

Thus, if extra dimensions do exist, and the fundamental scale of those dimensions is 
below the reheat temperature, a relic thermal CGWB ought to exist today. Compared 
to the relic thermal photon background (CMB), a thermal CGWB would have the same 
shape, statistics, and high degree of isotropy and homogeneity. The energy density (p g ) 
and fractional energy density (Q g ) of a thermal CGWB are 



,r 2 /-3.91^ 4/S 



Y 



15 V 9* ) 



Tcmbi (46) 



n = El ~ 3.1 x io-%0 4 / 3 , (47) 

Pc 

where p c is the critical energy density today, Tomb is the present temperature of the 
CMB, and g* is the number of relativistic degrees of freedom at the scale of mo. g* 
is dependent on the particle content of the universe, i.e. whether (and at what scale) 
the universe is supersymmetric, has a KK tower, etc. Other quantities, such as the 
temperature (T), peak frequency (v), number density (n), and entropy density (s) of 
the thermal CGWB can be derived from the CMB if g+ is known, as 

/3.91\ /3.91\ 
n g = n C MB f — J , s g = scmb f — J (48) 

3.91\ 1/3 /3.91 x 1/3 



T ' = TcMB \—j ' = VCMB \17) ' (49) 

These quantities are not dependent on the number of extra dimensions, as the large 
discrepancy in size between the three large spatial dimensions and the D extra 
dimensions suppresses those corrections by at least a factor of ~ 10~ 29 . If mo is just 
barely above the scale of the standard model, then g*= 106.75. The thermal CGWB 
then has a temperature of 0.905 Kelvin, a peak frequency of 19 GHz [I]. 
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